Dynamics of a viscous vesicle in linear flows 
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An analytical theory is developed to describe the dynamics of a closed lipid bilayer membrane 
(vesicle) freely suspended in a general linear flow. Considering a nearly spherical shape, the solution 
to the creeping-flow equations is obtained as a regular perturbation expansion in the excess area. 
The analysis takes into account the membrane fluidity, incompressibility and resistance to bending. 
The constraint for a fixed total area leads to a non-linear shape evolution equation at leading order. 
As a result two regimes of vesicle behavior, tank-treading and tumbling, are predicted depending on 
the viscosity contrast between interior and exterior fluid. Below a critical viscosity contrast, which 
depends on the excess area, the vesicle deforms into a tank-treading ellipsoid, whose orientation 
angle with respect to the flow direction is independent of the membrane bending rigidity. In the 
tumbling regime, the vesicle exhibits periodic shape deformations with a frequency that increases 
with the viscosity contrast. Non-Newtonian rheology such as normal stresses is predicted for a dilute 
suspension of vesicles. The theory is in good agreement with published experimental data for vesicle 
behavior in simple shear flow. 
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PACS numbers: 47.15.G-, 83.50.-v, 87.16Dg 

I. INTRODUCTION 

The dynamics of deformable particles such as drops, 
capsules and cells in flow represents a long-standing 
problem of interest in many branches of science and en- 
gineering, for instance because of its relevance to the 
rheology of emulsions and biological suspensions such as 
blood. The problem is challenging because the shape of 
these "soft" particles is not given a priori but is gov- 
erned by the balance between interfacial forces, e.g. due 
to stretching and/or bending of the interface, and fluid 
stresses. The interfacial properties, therefore, play a cru- 
cial role in the dynamics of these particles. 

The interface between two simple fluids is governed by 
the surface tension, which is isotropic, in the absence of 
surfactants or heating. Surface tension acts to minimize 
the interfacial area; therefore the rest shape of a drop is a 
sphere. Under shear flow, the drop deforms initially into 
an ellipsoid. As the flow strength increases, the drop area 
also increases and the angle between the ellipsoid major 
axis and flow direction decreases from 7r/4. If the flow 
strength is sufficiently high, and the viscosity contrast is 
moderate, the drop breaks up. Drop microhydrodynam- 
ics has been studied quite extensively^. 

Polymerized membrane interfaces found in synthetic 
capsules and the red blood cell, where a lipid bilayer 
is attached to a scaffolding of cytoskeletal proteins, ex- 
hibit far more complicated mechanical properties. These 
membranes behave as thin viscoelastic materials that 
can develop bending moments similar to thin elastic 
shells. Barthes-Biesel^ has reviewed the various consti- 
tutive laws adopted to describe the membrane mechanics 
and the effects of interfacial properties on the capsule be- 
havior in a linear flow. Yet, most of the theoretical stud- 
ies treat the membrane as a two-dimensional viscoelastic 
surface with no bending resistance. 

Bending stiffness, however, plays a crucial role in the 
mechanics of biological membranes, for example the equi- 



librium biconcave shape of the red blood cell can not 
be accounted for without including bending resistance^. 
The main structural component of the cell membrane is 
the lipid bilayer and its mechanical properties are es- 
sential to the overall cell mechanics. The pure lipid 
membrane consists of two sheets of lipid molecules. The 
molecular thickness imparts bending resistance. Lipid 
molecules are free to move within the monolayer, and 
therefore, in contrast to solid-like polymerized mem- 
branes, the lipid bilayer membrane is fluid. In addi- 
tion, since the lipid bilayer contains a fixed number of 
molecules with fixed area per molecule (under moderate 
stresses), the membrane is incompressible and the total 
area is constant. The enclosed volume is also constant at 
given osmotic conditions. The mechanics of lipid bilayers 
in concisely reviewed in Powers^. 

Equilibrium mechanical properties of vesicles made of 
lipid bilayer membranes are fairly well understood 6 , and 
vesicle shapes can be generated by minimizing the mem- 
brane bending energy subject to the constraints of con- 
stant area and enclosed volume. In contrast, the non- 
equilibrium dynamics of lipid bilayer membranes has 
been studied only to a limited extent. Experimental 
studies£&2JiI of vesicle behavior in unbounded shear flow 
observe that in weak flows, and when the inner and outer 
fluids are the same, the vesicle deforms into a tank- 
treading stationary prolate ellipsoid with an inclination 
angle close to 7r/4 with respect to the flow direction; 
however, in striking contrast to drops, when the fluid 
inside is more viscous than outside the vesicle under- 
goes a tumbling motion. Numerical simulat ion o 11 ' 12 ! 13 
and analytical theorie o 14 ' 15 ! 16 of vesicle microhydrody- 
namics attempt to elucidate such experimental observa- 
tions. In their classic paper Keller and Skalakii analyzed 
the motion of a tank-treading ellipsoid in unbounded 
shear flow. The theory qualitatively captures features 
like the tumbling transition, but their assumption of a 
fixed particle shape casts some doubts on the applica- 



2 



bility to deformable vesicles; for instance, no connection 
can be made to the physical properties of real membranes 
such as the bending rigidity. The free-boundary charac- 
ter of the problem is taken into account in several recent 
works. Shape evolution of a fluctuating quasi-spherical 
vesicle was considered by Seifert^ for the case where 
the inner and suspending fluids are the same, i.e. there 
is no viscosity contrast. The theory predicts a station- 
ary, tank-treading prolate ellipsoid and no transition to 
tumbling motion. The importance of viscosity contrast 
for the tumbling transition was recognized by Misbah 16 , 
who showed quantitatively that a critical value of viscos- 
ity contrast exists which separates the tank-treading and 
tumbling regimes. The critical viscosity ratio was shown 
to decrease with the excess area (the difference in the ar- 
eas of the vesicle and an equivalent sphere with the same 
volume). He also pointed out that the area constraint 
leads to non-linear leading order evolution equations for 
the vesicle shape, which in the case when only ellipsoidal 
deformation modes are considered are independent of the 
membrane bending rigidity. Earlier work by Olla 15 de- 
rived similar results for a viscoelastic membrane. 

Cell behavior in unbounded flow is of fundamental in- 
terest. However, flows in confined geometries are much 
more relevant to physiology, for instance blood flow in 
the microcirculation. Walls affect greatly particle mi- 
crohydrodynamics, for example red blood cells migrate 
away from the blood vessel walls, an observation that 
dates back to Poiseuille 1 -. The existence of a near-wall 
cell-depleted region accounts for the Fahraeus-Lindqvist 
effect, which is the decrease of the apparent blood vis- 
cosity in smaller vessels^. Cell traffic between the blood 
stream and tissues involves cell attachment to blood 
vessel walls^; examples are leukocytes during inflam- 
matory response, platelets in formation of atheroscle- 
rotic plaques, or tumor cells in metastasis. In order to 
elucidate the fundamental features of the process, sev- 
eral studies have used vesicles as the simplest artificial 
cel l 21 ' 22 ! 23 . They have reported that the flow-induced 
deformation of adhering vesicles gives rise to a lift force 
that can lead to unbinding from the substrate. The ex- 
perimental data on the lift force and vesicle migration 
velocity have not yet been quantitatively compared to 
theoretical studies 2 ^ 5 .^. 

The purpose of this paper is to simultaneously include 
the effects of (i) viscosity contrast, (ii) membrane in- 
compressibility and (Hi) bending rigidity in the analysis 
of vesicle dynamics in shear flows, unbounded or in the 
presence of a wall, in a consistent way and thereby pro- 
ceed further towards a fully quantitative description of 
the experimentally observed vesicle behavior. As a first 
step, the leading order small deformation of a nearly- 
spherical vesicle will be considered. However, the devel- 
oped formalism will serve as a rigorous basis for consider- 
ing higher orders in the non-linear dynamics of a vesicle 
in flow. 

Our study extends the results of Seiferti^ derived for an 
equiviscous vesicle. The works by Misbahi^ and Olla 15 



are clarified, more specifically, the calculations are pre- 
sented more explicitly, some important results are cor- 
rected and cast in a form that can be directly compared 
to experimental measurements, and more physical in- 
sights are given for the vesicle behavior in the tumbling 
regime. We demonstrate that the theory agrees quantita- 
tively with published experimental dat a 8 ' 9 ' 10 . We present 
new results for (i) the non-Newtonian rheology of a vesi- 
cle suspension, in particular, we show how the single vesi- 
cle solution serves to calculate the effective stress of a 
collection of many vesicles and we predict the existence 
of normal stresses, and (ii) vesicle migration velocity in 
wall-bounded flow in the case of large distances from the 
boundary. 



II. PROBLEM FORMULATION 

Let us consider a neutrally-buoyant vesicle formed by 
a closed lipid bilayer membrane with total area A. The 
vesicle is suspended in a fluid of viscosity 77 and filled with 
a fluid of viscosity (A — 1) r\ . Both interior and exterior 
fluids are incompressible and Newtonian. The vesicle has 
a characteristic size a defined by the radius of a sphere 
of the same volume. The equilibrium shape of the vesicle 
is characterized by a small excess area 

A = A/a 2 -4:TT. (1) 

The coordinate system employed is spherical (r, 9, 4>), 
with the origin coinciding with the center of mass of the 
vesicle. The interface is specified by a shape function 
F(r,t), of the position r and time t, which represents 
the interface as the set of points, where F (r, t) = 0. It 
is given by the relation 

F(r,t)=r-r e (0,4>,t). (2) 



A. Hydrodynamics 

The vesicle is placed in a steady two-dimensional linear 
flow 

v°°(r)= 7 E-r, (3) 

where 7 is the strain rate, and E is the velocity gradient 
tensor. Linear flows are defined by 

, / 1 + w 0\ 
E=- [l-u , (4) 
1 \ / 

where ui is the magnitude of the rotational flow com- 
ponent. Irrotational flow such as the pure extensional 
flow is given by oj — 0. Simple shear flow is specified by 
u> = 1, i.e. =72/. A sketch of the problem is shown 
in Figure [TJ 
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where k is the bending rigidity, H and K are the mean 
and Gaussian curvatures, respectively, and a is the local 
membrane tension. The tangential part of the surface 
force density arises from nonuniformitics in the surface 
tension, which are needed to ensure local area incom- 
pressibility, as discussed in more details in the next sec- 
tion. The mean and Gaussian curvatures are 



(11) 



FIG. 1: A sketch of a vesicle is a simple shear flow. The zoomed 
region of the interface illustrates the bilayer lipid structure of the 
membrane. 



Typically vesicles are micron-sized. At these small 
length-scales water is effectively very viscous, inertial ef- 
fects are unimportant and fluid velocity fields inside v ln 
and outside v out the vesicle are described by the Stokes 
equations 



7yV 2 v out 



0, v. 



0, 



(A - l^V V n - V P iu = 0, V • v in = . 



(5) 



Far away from the vesicle, the flow field tends to the 
unperturbed external flow v out — > v°°. Velocity is con- 
tinuous across the interface 



at 



(6) 



where r s denotes the position of the interface. The inter- 
face moves with the fluid velocity-21 



dF 
~dt 



v s • VF = . 



(7) 



where v s is the fluid velocity at the vesicle interface. 
Fluid motion gives rise to bulk hydrodynamic stress 



-pi + r)(Vu + (Vu) 



(8) 



where I denotes the unit tensor and the superscript T 
denotes transpose. The surrounding fluids exert tractions 
on the membrane that are balanced by membrane forces 



n 



at 



(9) 



where n is the outward unit normal vector and the mem- 
brane surface forces are discussed next. 



B. Membrane mechanics 

Unlike drops, which are governed by surface ten- 
sion, the shape of the vesicle is determined by bending 
elasticity 6 . A nonequilibrium membrane configuration 
gives rise to a surface force density^ 



[2aH - k {AH 3 
- V s cr , 



4KH + 2W„H) \ n 



(10) 



K = -V • (nVn + n x (V x n)) , (12) 

where the outward unit normal vector can be determined 
from the shape function ([2]) via 



VF_ 

Wf\ 



(13) 



The surface gradient operator is defined as 

V S = I S -V, (14) 

where the matrix I s = I — nn represents a surface projec- 
tion. At rest, Eq. ([9]) reduces to the well-known Euler- 
Lagrange equation 

p in _ p out = 2(jH _ K \^ H z _ iKH + 2\/ 2 s H] . (15) 

C. Area constraint and tension 

The lipid bilayer is fluid because the lipids can dif- 
fuse rapidly within the membrane. Moreover, since the 
number of lipids in a monolayer and the area per lipid 
are fixed, the lipid bilayer membrane is incompressible 
and the total area is conserved. A membrane element 
only deforms but does not change its area. Accordingly, 
the local tension changes in order to keep the local area 
constant. Hence, inhomogeneities in the tension ensure 
local area incompressibility. The situation is analogous 
to three-dimensional incompressible fluids, where finite 
changes in pressure (tension is the two-dimensional ana- 
logue) correspond to infinitesimal changes in fluid den- 
sity, and pressure takes the place of the density as an 
independent field variable; in flowing fluids pressure can 
become nonuniform. The local area conservation implies 
that the velocity field at the interface is solenoidal 14 



V, • v = . 



(16) 



The global area constraint acts like an isotropic ten- 
sion whose value at equilibrium is given by the Lagrange 
multiplier used to determine the shaped. Under flow, 
the changes in shape result in variations of this effective 
isotropic tension. 
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D. Dimensionless parameters 

Viscous forces exerted by the extensional component 
of the external flow drive shape deformation that occurs 
on a time scale 



A7 



(17) 



Several intrinsic relaxation mechanisms oppose the defor- 
mation. Bending stresses work to bring the shape back 
to its preferred curvature state; the corresponding time 
scale is 



t K 



Xrja 3 



(18) 



In shear flow, vesicle rotation away from the extensional 
axis of the imposed flow effectively decreases the extent 
of shape distortion; the associated time scale is 



(19) 



The strength of these two relaxation mechanisms that 
limit shape deformation by the flow is quantified by the 
corresponding dimensionless parameters: the capillary 
number 



x = — 



and the rotation parameter 



= \- 



(20) 



(21) 



The smaller of these parameters controls the magni- 
tude of vesicle deformation. At moderate viscosity ra- 
tios, the vesicle shape should remain close to the equi- 
librium one provided that the capillary number is small, 
i.e. restoring bending forces are stronger than the dis- 
torting viscous forces. A high viscosity inner fluid limits 
the shape distortion in flows with non-zero vorticity by 
means of increasing the rate of vesicle rotation. In addi- 
tion to the flow-related parameters, a physical parameter 
that arises from vesicle geometry is the excess area ((T|), 
which sets the maximum magnitude of the shape distor- 
tion; it turns out to be the relevant small parameter for 
the analysis of nearly-spherical vesicles as shown in the 
next sections. 

Henceforth, bending stresses and tension are normal- 
ized by k/o 2 ; all other quantities are rescaled using i], a, 
and 7. Accordingly, the time scale is 7 -1 , the velocity 
scale is 7a, bulk stresses are scaled with 777. 



III. SMALL DEFORMATION THEORY 

In this section we present a perturbative solution for 
the microhydrodynamics of a vesicle with a nearly spher- 
ical shape. 



In the coordinate system centered at the vesicle, the ra- 
dial position r s of the vesicle interface can be represented 
as 



r s = a + f (SI) , 



(22) 



where / is the deviation of vesicle shape from a sphere, 
which depends only on the angles (9 , 4>) (or equivalently 
the solid angle SI) and has a vanishing angular average 



fdSl = 0. 



(23) 



The isotropic contribution a is determined by the 
volume-conservation constraint 



J (a + f) 3 dSl = 4tt . 



(24) 



The total area conservation constraint relates the ampli- 
tude of the perturbation / and the excess area A 



-dSl = 47r + A . 



A/a 2 = 



where f denotes the unit radial vector. 



(25) 



A. Expansion in spherical harmonics 

In Eq. (j2"2")l , the function / representing the perturba- 
tions of the vesicle shape depends only on angular coordi- 
nates. Thus, it is expanded into series of scalar spherical 
harmonics Yj m (IA1|1 



00 j 



f — fjmYjm , 

j=2 m=-j 



(26) 



In the above equation, the summation starts from 
nonzero j because / includes only the nonisotropic contri- 
butions. For small shape perturbations around a sphere 
the volume constraint (l24l) become a 14 ' 29 



(27) 



jm 



where the sum over j starts from 2, \m\ < j and 
f* m = (— l) m fj_ m . Similarly, the area constraint 
transforms to 



where 



U + 2) (j - 1) 



(28a) 



(28b) 



The tension is expanded in scalar spherical harmonics 

00 j 



j=2 m=-j 



(29) 
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where <7o is the isotropic part of the tension, which varies 
with shape in order to keep the total area constant. It 
is determined from the condition that the modes must 
satisfy the area constraint A = 0, i.e. 

E a ^™/;™=°- ( 3 °) 

jm 

The nonuniform part of the tension is related to the local 
incomprcssibility. 

B. Perturbation solution 

The combination of perturbative analysis and the 
spherical harmonics formalism for solving Stokes-flow 
problems involving a deformable particle has been devel- 
oped in detail in Vlahovska et aZi 29 ' 30 for the problem of 
a surfactant-covered drop. Here we present a brief out- 
line of the method, and focus on the new feature specific 
to lipid bilayer membranes which is the area constraint. 

This paper considers a non-spherical particle albeit the 
shape deviation from sphericity is small, i.e. / ~ 0(e) 
and e<1, In this case, the exact position of the inter- 
face is replaced by the surface of a sphere of equivalent 
volume, and all quantities that are to be evaluated at 
the interface of the deformed particle are approximated 
using a Taylor series expansion. 

1. Small parameter for area-constrained dynamics 

Small deformations of initially spherical drops and cap- 
sules with elastic membranes have been considered in a 



f jm = iw — f jm + A 1 C(X,j 1 m)+\ 



where the isotropic tension cro is determined using the 
area constraint (|28p as described in Appendix IFl see Eq. 
(|F4[) . The first term in the evolution equation (|32|) de- 
scribes rigid body rotation of a particle with shape / 
by the rotational component of the external flow. The 
second term describes the distortion of the vesicle shape 
by the cxtcnsional component of the external flow. The 
term including x is associated with relaxation driven by 
the membrane stresses. The expressions for C and T are 
given by (|E12|) . These coefficients depend on A and are 
bounded at A — > oo. 

The area constraint couples all modes, and results 
in isotropic tension, which depends nonlinearly on the 
shape. Thus the leading order vesicle dynamics is non- 
linear in contrast to the corresponding results for drops 
and capsules. As we show in the next section, a peculiar 



number of studie s 31 ' 32 . Typical choices for the small pa- 
rameter are the capillary number or the inverse viscosity 
ratio. 

In contrast to drops and capsules, the rest shape of 
the vesicle we consider is non-spherical. The excess area 
plays a crucial role in the vesicle dynamics, because under 
the constraints of constant area and volume a sphere is a 
geometrically rigid object. Consequently, in the absence 
of excess area, a vesicle behaves as a rigid sphere. 

Thus, the appropriate small parameter that reflects the 
importance of the excess area in the vesicle dynamics is 

e = A^. (31) 

The square root comes from the observation that f 2 ~ 
A (|28p. The choice of the excess area as the small pa- 
rameter allows vesicle behavior to described, where the 
whole excess area is involved in the vesicle deformation. 



2. Evolution equation 

We solve the hydrodynamic problem to obtain the ve- 
locity field, and use the fact that the interface moves with 
the normal component of the velocity to determine the 
shape evolution. The solution for small deviations from 
a sphere is presented in details in Appendix [EJ In this 
section we present the general expression for the shape 
evolution equation valid for any linear flow. In the subse- 
quent sections we analyze in more details the particular 
case of simple shear flow. 

At leading order the evolution of the shape parameters 
(f2"o) is described by 



1 X- 1 r(A,ao,j)/ jm + 0(A- 1 £,£ 2 ) , (32) 
I 

consequence from the area-constrained dynamics is that 
the stationary solution is independent of the capillary 
number. 



IV. RESULTS: SIMPLE SHEAR FLOW 

In this section we provide analytical solutions for the 
vesicle shape evolution in a simple shear flow. Sim- 
ple shear flow consists of an extensional component, 
v ext ~ \{Vi x ify : which in the spherical harmonics rep- 
resentation (|D1|) is described by 

c 2±20 = ~Fi\/f c 2±22 = "Fiy ff , (33a) 
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and a rotational component, vj^ t 
specified by 



oo 
-101 



h(y, —x, 0), which is 



(33b) 



The extensional part of external flow ([3]), which is re- 
sponsible for shape distortion, is fully characterized by 
a second-order traceless tensor, which corresponds to 
spherical harmonics of the order j = 2 (|33a[) . Therefore, 
at leading order the shear flow affects only the subspace 
j = 2. Considering only these modes simplifies the ex- 
pression for the tension (|F4[) and the evolution equations 
(13"!?)) to the following set of coupled non-linear differential 
equations 



h 



2 



2ifo(A)A- 1 (/ 22 -/ 2 _ 2 )/ 2 _, 



(34) 




A 1/2 (9 + 23A)/(16tt 3/2 > 



FIG. 2: Stationary inclination angle in simple shear flow as func- 
tion of the rescaled excess area for different viscosity contrasts. The 
symbols are experimental data from Kantsler and Steinberg^. The 
lines represent the linear small deformation theory (I4U1 . 



where 



h(X) 



23A + 9 



(35) 



Strictly speaking Eq. (|3~4"|) is only valid at long times 
(t ^> t K ) when all transient j ^ 2 modes have decayed. 
The initial conditions are set by the equilibrium vesi- 
cle configuration. For example, for a fluctuating quasi- 
spherical vesicle^ 



1/ 



jm\ 



KE(j, a eq ) ' 



(36) 



where E(j, a) is given by (|E12d|) . 

The shape parameters (f2l)]) are decomposed into real 
and imaginary parts 



f ■ = f'. 4- i f " f* = f. — i f" 

Jim Jim ' Jim Jim Jim J nr. 



(37) 



In the flow plane x—y, the vesicle shape / is characterized 
by only three components, / 22 ,/ 22 , and / 2 o, correspond- 
ing to deformation along the flow axis x, straining axis 
x = y, and the z axis. Shape modes / 2 ±i describe defor- 
mations of the type xz and yz, which vanish in the flow 
plane z — 0. From geometrical considerations we have 
for the inclination angle 



b = -i arctan ( 

1 \J22 



(38) 



The numerical solution of the full set of evolution equa- 
tions (|3"2"|) shows that all other modes decay to zero. 
Hence, in the stationary state all excess area is stored 
in the / 2 ± 2 modes. Substituting the stationary ampli- 
tudes (|39| in the relation for the inclination angle (|38)) 
and expanding for small values of the excess area A we 
obtain 



J o = ~r 



(9 + 23A)A 



1/2 



(40) 



which agrees with the expression reported by Seiferti^ 
for the case of no viscosity contrast (A = 2). The analo- 
gous result reported by Misbahi^ contains a misprint (140 
should read 240). Our relation (j40|) is in good quantita- 
tive agreement with the experimental data of Kantsler 
and Steinberg 8 , as demonstrated in Figure [2] Eq. (f39|) 



implies that if A > Ah , i.e. if the viscosity contrast is 
higher than a critical value 



9 120 

c ~~ ~23 



2tt 



23 V 15A 



(41) 



no stationary tank-treading solution exists, as also 
pointed out by Misba h 16 i 33 , and the vesicle starts to tum- 
ble. 



B. Tumbling 



A. Tank— treading 



The evolution equations (|34|) have two sets of sta- 
tionary points. The first one corresponds to the tank- 
treading state. The only non-zero stationary amplitudes 
are 



■til 



4ft(A) 



Ah 2 (\) 



4/i(A) 



(39) 



The second fixed point of Eq. (f3~4")l corresponds to peri- 
odic vesicle deformation given by oscillating j = 2 modes 
around values 



f22=KX), #2=0 /^=0, /»=0. 



(42) 



If /2o(0) 7^ 0, the J 2 o(i) is also oscillating. All modes 
with j ^ 2 decay to zero and are not oscillating. 

The time-periodic vesicle deformation depends 
strongly on the viscosity contrast as illustrated in 
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of the major and minor axes 



FIG. 3: Time dependence of the /2 m modes in the (a) breathing 
(A > A c ) and (b) tumbling (A S> A c ) regimes; solid, dashed and 
dotted lines represent the f^, an d /20 modes, respectively. The 
excess area is A = 0.2. The viscosity contrasts are (a) A = 10, (b) 
A = 100. The initial conditions are / 22 (0) = /22(D) = VO.IA; 
/20(D) is determined by the area constraint. 



Figure [3j Figure [3] (a) shows that close to the critical 
viscosity contrast the period of oscillations is long and 
the amplitude of oscillations of the /20 mode, which 
corresponds to out-of-the-flow-plane deformation, is 
significant. Experimentally the vesicle appears to be 
trembling or "breathing" 9,10 . The "breathing" motion 
has been discussed in parts by Misbahi^ although the 
role of the /20 mode has not been recognized. 

The mode oscillation frequency increases with viscos- 
ity contrast, as seen in Figure [3] (b). At high viscosity 
contrast it approaches twice the rate of rotation of the 
external flow, 7. This can be seen from Eq. (f3~4"| where 
in the limit A ^> 1 only the rotation terms survive. Since 
the /20 mode is not rotationally stabilized, its oscillation 
amplitude becomes negligible at high viscosity contrast; 
correspondingly the out-of-the-flow-plane deformation is 
suppressed. The periodic vesicle deformation then corre- 
sponds to rigid body rotation, /^(t) ~ cos(t) , W ~ 
sin(i). 

In order to compare theory and experiment we derive 
expressions for some experimentally measurable param- 
eters. Experimentally only the deformation in the flow 
plane is observable. The vesicle contour represents an el- 
lipse. Combining (J^U) and (|38|) we obtain for the lengths 



' max,rmn 



(t) = 1 - f W (t) 



± A 



ii 

2tt ' 



where we have neglected the contributions from the /2±i 
modes for the sake of simplicity. We see that the major 
and minor axes follow the oscillations of the /20 mode: 
the larger the /20 oscillation amplitude, the larger the 
amplitude in the r max and r m i n fluctuations. In the ex- 
periments of Kantsler and Steinberg 9 -, the excess area A 
corresponding to a prolate ellipsoidal shape with major 
and minor axes as seen in the flow plane was reported. 
Using the relation between the deformation in the flow 



plane and the excess area (r ma x - r m in)l {r ma x + r„ 
(15A/327r)2, reported first by Seifertii, we obtain 



A(i) = 



A - 2fj Q (t) 
1 - ho{t)\[^ 







(44) 



Unlike the true excess area A, which is constant, A is 
time-dependent due to the fluctuations in the /20 mode, 
i.e. the out-of-the-flow-plane transfer of area. 
The angle evolution is described by 



h(X) 



(A 



3(20o (t)), (45) 



from which we infer that the tumbling period, defined 
as the time needed for a material point to return to its 
initial position, is given by 



Ti 



tumble 



= in 



1 - 



4/i(A) 5 



(46) 



As already discussed, when the viscosity contrast is high 
the oscillations of the /20 mode are suppressed and 
/2o(i) ~ /2o(0). In this case Eq. (|4"S"j) becomes equiva- 
lent to the Keller-Skalak equation <fio = A + B cos(2</>o), 
where the constants A and B are identified with —1/2 
and h(X) / \J A — 2/| (0) , respectively. 

We compared the time-dependent vesicle behavior 
with the experiments of Mader et al^. A good agree- 
ment was obtained for the angle evolution, for example 
of their vesicle number 6 using reduced volume 0.988 in- 
stead of the reported 0.996. The difference is reason- 
able given the uncertainty with which the reduced vol- 
ume is determined experimentally: it is computed as- 
suming axial symmetry about the vesicle's longest axis, 
which might not be always the case because of the pres- 
ence of a non-zero /20 mode. Note, however, that in the 
tank-treading regime, since the excess area is stored only 
in the /2±2 modes the vesicle is a prolate ellipsoid and 
the computation of the reduced volume is quite accurate. 
We also attempted to reproduce theoretically the exper- 
imental curves in Figure 4 of Kantsler and Steinberg 9 . 
The predicted tumbling frequency from Eq. (|46p using 
viscosity contrast and excess area close to the reported 
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values was almost twice those of the experimental data. 
One possible explanation is that our small-deformation 
theory fails for the large values of excess area in this ex- 
periment (A« 1). 



C. Rheology 

A suspension of vesicles can be described as a contin- 
uum with effective properties at length scales large com- 
pared to the size of the constituent particles. In dilute 
suspensions particles are far from each other and they do 
not feel each others presence. Thus, the effective stress 
is just a sum of the stress contributions of the individ- 
ual particles and the bulk stress of a dilute suspension is 
linear with the particle volume fraction y i 34 i 35 



£ = 2W + <pT , 



(47) 



where E s is the symmetric part of the velocity gradient 
tensor ([4]) and T is the stress associated with the vesicle 
disturbance velocity field. Rheological properties of in- 
terest are the particle contribution to the shear viscosity, 
T xy , and the normal stress differences, N\ — T xx — T yy 
and N 2 = T yy — T zz . Vesicle stress is directly related 
to the amplitudes of the velocity field j = 2 (jDlji — cor- 
responding to the stresslet (the symmetric force dipole). 
Using the relations (|C1[) we obtain for the vesicle contri- 
bution to the effective shear viscosity 



T 



5 23A 
2 23A+9 



39 



X 



L 3i 



30 £(2,<Tq)(/22-/2-2) 

23A+9 



(48) 



where the shape deformation modes fa±2 are given by 
(|E13|) . For a given excess area and viscosity contrast the 
effective viscosity in the tank-treading regime is 



T^_ A (23A + 9) 
v 2 16tt 



(49) 



which is obtained by inserting the expression for the 
modes amplitude (|E13|) in Eq. (148|) and taking into ac- 
count (|E15|) . In the limit of a spherical particle, A = 0, 
and since a sphere with fixed area and volume in shear 
flow undergoes only a rigid body rotation, the Einstein's 
result for a suspension of rigid spheres is recovered. The 
effective viscosity decreases with the increase of the ex- 
cess area because the deformable vesicles elongate and 
thus offer less resistance to the flow. An increasing vis- 
cosity contrast also leads to a decrease in the effective 
viscosity because vesicles align better with the flow (the 
inclination angle (I40[) decreases with viscosity contrast). 

Unlike a suspension of rigid spheres, a suspension of 
deformable vesicles exhibits normal stresses. In the tank- 
treading regime below the critical viscosity contrast, the 
magnitudes of the first normal stress difference is 



Ni = 



A(23A + 9) 
8^ 



-1 



1920tt 



A(23A + 9) S 



(50) 



D. Vesicle migration in the presence of a wall 

A spherical particle in simple shear flow produces a 
symmetric disturbance velocity field and therefore the 
particle does not drift relative to a bounding wall 37 . Par- 
ticle deformation breaks the symmetry and may lead to 
cross-stream migration. 

The leading order term in the far field of the distur- 
bance velocity due to a force-free and torque-free particle 
is the stresslet, the symmetric and traceless force dipole. 
The boundary conditions at the wall can be satisfied by 
placing the stresslet hydrodynamic image on the oppo- 
site side 3 ^. Thus, a particle far away from the wall moves 
with a velocity due to its corresponding image stresslet; 
in particular, the vesicle drift velocity normal to a rigid 
wall is proportional to the stresslet component in the di- 
rection of the plane unit normal, i.e. T yy . Smart and 
Leighton— report that 



16di Tytn 



(51) 



where d denotes the distance from the particle center 
to the wall. Taking into account that the normal stress 
is T yy — — (JVi — AT 2 )/3 and using the normal stresses 
expressions (T5TJ|) we obtain for the lift velocity 



where 



c(A,A) 



c(A,A) 
d 2 



32 



(52) 



(53) 



and the second normal difference is N 2 



Similar type expressions have been reported in analyti- 
cal studie d 26 ! 39 , although the dependence on the excess 
area has not been presented explicitly. Eq. (|53|) shows 
that the more deformable the vesicle, i.e. the larger the 
excess area, the larger the migration velocity. Interest- 
ingly, at very small excess area we obtain from ({50")) that 
Nt ~ Aa, and vesicle migration is independent of the 
viscosity contrast. The values of the prefactor ([55)) are 
of the same order of magnitude as reported in studies 
considering a tank-treading ellipsoid^ 9 - or numerical sim- 
ulations of a vesicle with no viscosity contrast. 



V. SUMMARY AND OUTLOOK 

This study considers the behavior of a vesicle formed 
by a fluid membrane in a general shear flow. The analysis 
takes into account the membrane incompressibility and 
bending resistance as well as the viscosity contrast be- 
tween the interior and exterior fluids. Analytical results 
for a nearly spherical vesicle are obtained. Expressions 
for the shape evolution equation, stress and velocity fields 
and the effective membrane tension, which is needed to 
enforce the area constraint, were derived using a formal- 
ism based on spherical harmonics. The main theoretical 
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result of this study is contained in the shape evolution 
equation Eq. (f3"2")l and the expression for the effective ten- 
sion Eq. (|F4|) . The latter serves to eliminate the tension 
in favor of the excess area, which is the physically relevant 
parameter. In contrast to drops and capsules, the shape 
evolution of the area-constrained vesicle is non-linear. 

The derived shape evolution equation describes vesicle 
dynamics in general linear flow. In the particular case of 
simple shear flow, only the shape modes with the same 
spherical harmonic order as the external flow, j ' = 2, are 
perturbed; all other modes decay on a time scale set by 
the bending rigidity. Consequently, the shape evolution 
equations simplify to Eq. ([34]) and yield expressions for 
quantities that can be measured experimentally such as 
the inclination angle with respect to the flow direction, 
the lengths of the major and minor axes of the vesicle 
contour in the flow plane, etc. In the stationary tank- 
treading state, we show that in a simple shear flow the 
leading order deformation and stresses are independent 
of the membrane bending rigidity. Our theory is in quan- 
titative agreement with experimental data for the vesicle 
deformation in the tank-treading and in the tumbling 
regimes. Non-Newtonian rheology with normal stresses 
is predicted for a suspension of vesicles. We also consid- 
ered the vesicle dynamics in wall-bounded shear flow and 
presented a simple derivation of the leading order correc- 
tion, in the distance to the wall, to the rate of vesicle 
migration away from the wall. 

Several problems remain open. We discussed the main 
features of the time-dependent vesicle dynamics in lin- 
ear flows, but this problem remains to be systemati- 
cally explored. Although the quantitative agreement be- 
tween predicted and experimentally observed 1 ^ behavior 
of tumbling vesicles with small excess area is encourag- 
ingly good, for vesicles with large excess area£ the agree- 
ment is poor, which might be due to the limitations of 
the leading order theory. Our analysis will be extended 
to consider higher-order perturbations, and hence eluci- 
date the feedback from the flow on vesicle dynamics. 

The developed theory applies only to a nearly spherical 
vesicle. Larger deformations from sphericity can be ex- 
plored by numerical simulations. A body of work exists 
on capsule dynamics, mostly considering elastic mem- 
branes with no bending resistanc o 40 ' 41 ! 42 . To our knowl- 
edge there are only a couple of numerical studies that in- 
clude bending resistance 43,44 . Area incompressible fluid 
membranes with bending resistance such as the lipid bi- 
layer membranes, have been studied only to a limited 
exten t 11 ! 12 ! 45 and some results are conflicting. For ex- 
ample, Noguchi and Gompper— report stationary tank- 
treading discocyte shape, while Kraus el alr^ find only 
prolate ones. More efficient and accurate simulations are 
needed in order to perform a systematic study on flow in- 
duced shape transitions. In order to explore numerically 
the dynamics of highly non-spherical vesicle shapes we 
are developing Boundary Integral Method simulations, 
which combine an algorithm for adaptive restructuring 
of the computational grid that allows resolution of high 



curvature regions 4 ^ with the area constraint (|16p and in- 
tcrfacial force density (TTUJ). 
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APPENDIX A: SPHERICAL HARMONICS 

For the sake of completeness, we list the definitions of 
scalar and vector spherical harmonic a 47 ! 48 . The normal- 
ized spherical scalar harmonics are defined as 



■^j'm (^) 



2j+l (j-m)\ 
4tt (J +m) ! 



(-l) m P™(cos6»)e imv , (Al) 



where f = r/r, [r,6, (p) are the spherical coordinates, 
and P.™ (cos 0) are the Legendre polynomials. The vector 
spherical harmonics are defined as 



Yjmo = [j U + 1)] 5 rV n Y jm , 



Yjml 



YjmO ■ 



Yjm2 



rY 



(A2a) 



(A2b) 



(A2c) 



where Vn denotes the angular part of the gradient oper- 
ator, yjmo and y jm i are tangential, while y jm 2 is normal 
to a sphere. 



APPENDIX B: FUNDAMENTAL SET OF 
VELOCITY FIELDS 

Following the definitions given in Blawzdziewicz et 



i 36 



, we list the expressions for the functions u ■ (r, fi) 



u jmO = \ r 3 ( 2 - J + 3 r 2 ) YimO 

+fr-^(i + l)] 1/2 (l-r- 2 )y 3 > 



jml 



r 3 1 yjmi ■ 



u iro2 = ¥ 3 ( 2 - 3) (ifj) 1 (1 - r 2 ) Vimo (Blc) 
+ ir^(j + (2-j>- 2 )y, m2 . 



(Bla) 



(Bib) 
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< m o = ¥^ (-U + 1) + U + 3)r 2 ) y jm o , R9 . 
-i» J -i[,-( J - + l)] l/a (l-r>)y Jma , ^ 



(B2b) 



u jm2 ~~ 2' 



(3+i)(^) 1/2 (l-r 2 )y iro0 (B2c) 
+5^' _1 (i + 3 - (j + l)r 2 ) y jm2 ■ 



On a sphere r = 1 these velocity fields reduce to 

U f mq = yjmq • (B3) 

Hence, and u - ml are tangential, and u^ rl2 is normal 
to a sphere. In general, a vector velocity field which is 
tangential to a surface with normal n has an irrotational 
components* 

n • (V, x V rr ) = , (B4) 

and solenoidal component 

V s -v so '=0. (B5) 

On a sphere, the irrotational component is identified with 
the q = vector spherical harmonic, and the solenoidal 
corresponds to the q — 1 vector spherical harmonic. 



APPENDIX C: EFFECTIVE STRESS OF A 
DILUTE DISPERSION 

The disturbance velocity field due a particle can be 
represented as a superposition of velocity fields gener- 
ated by a collection of force multipoles, by analogy to 
electrostatic a 34 ' 50 . 

The strength of the stresslet, which is the symmetric 
and traceless force dipole moment, gives the particle con- 
tribution to the effective stress of a dilute dispersion 



T xy — — (C22O — C2-20) 



1 »r 3 / 1 . 

N 2 = --iVi + -\ — C200 ■ 

2 4 V 07T 



(Cla) 



(Clb) 



(Clc) 



The stresslet field is related to the amplitude of the j = 2 
velocity field; the relations between c and c are given by 



_ (2j+l)(2j-l)[j(j + l)] 1 / 2 

j(2j+l)(2j-l) 
T J -_ ) _ 1 t-jm2 ■ 



(C2) 



The complete expressions can be found in Blawzdziewicz 
et al^. 



APPENDIX D: VELOCITY FIELDS AND 
HYDRODYNAMICS STRESSES 

Velocity fields are described using basis sets of funda- 
mental solutions of the Stokes equations^, , defined 
in Appendix [Bl 

v° ut (r) [u+ m? (r)-uj mq (r)] 



jmq 



+ E/ C j<nq n jmq ( r ) ' 

v" 1 (r) = E C 3»'9 U « ( r ) 
jmij 



(Dla) 



(Dlb) 



The hydrodynamic tractions exerted on a surface with a 
normal vector n are n • T 



T = n T = TjmqYj 



mq 



(D2) 



In the particular case of a sphere characterized with a 
normal vector f , the viscous tractions are linearly related 
to the velocity field 



-out 
jmq 



E C TnW ( Q t'q - Q q'q) + E <W ^ (D3a) 



r in -Vc ,0+ 
q' 



(D3b) 



where 0*, are obtained from the velocity fields (|B1[) 



0') 



2j + i - 3 (i±iV\ 
^ i-i o J 

3 (^ii) 5 n •> ; 1 - : 



(D4) 



(i) 



2j + i + 4 y 

^ -2i-i 3(^y \ 

-j-2 



(D5) 

In the derivation of the above expressions we have used 
that 

f.(Vv + (Vv) r ) = r| (-Ulv(vr). (D6) 
dr \r J r 



APPENDIX E: LEADING ORDER 
PERTURBATION SOLUTION 

For small deviations from sphericity 2 ^, the 0(e) ex- 
pansions for the normal vector and the mean curvature 
are 



n (fi) = f - f jm [j {j + 1)] 2 y jm0 



(El) 
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h (n) = 1 + 1 (j + 2) (j - 1) / Jm r jm . (E2) 

At order 0(e), in the bending stresses (fTO|) the terms ff 3 
and iJX cancel and only V 2 _ff remains 



V 2 S H = - —j (j + 1) (j - 1) (j + 2) / Jm y jm . (E3) 

Combining Eqs. (|E1[) and (|E2|) in the Laplace term of 
the membrane stresses (TTTfl) yields 

oH = (i(7 (j - 1) (j + 2) / Jm + a jm ) Y jm . (E4) 

In the above equation the isotropic part has been omitted 
because it has no importance to flow dynamics; <Jj m fjm 
is neglected as well as it is a higher or der term. On a 
sphere V s cr becomes simply \fj{j + l)pjm~yjmO accord- 
ing to Eq. (|A2a|) . 

All modes that contribute to the excess area are af- 
fected by the flow. Therefore, we will present results for 
any j. In this way, the theory can be applied to fluc- 
tuating vesicles and other types of flow, e.g. quadratic 
(Poiseuille) flow. 

The incompressibility condition (Till)) implies that the 
amplitudes of the velocity disturbance field (|Dip are re- 
lated 



1 



Cjm2 = ^ VJ(J + !) c imO • 



(E5) 



At leading order the stress balance reads 

f • T out - (A - l)f • T in = t mcm . (E6) 



It gives a relation for the tractions (|D2[) . which is linear 
in j and m 



out _ (\ _ 1 \_in _ i 

jmq \" ) jmq 



mem 
jmq ■ 



(E7) 



The tangential interfacial stress has an "irrotational" 
component 



Linem „ — 1 



jmO 



jm ■ 



(E8) 



The normal component of the membrane stress is 



X- L (2a, m +E{j,a )f jm ) ■ (E9) 



irnem „ ,— 1 

z jm2 

The tangential balance (|E8j) determines the tension dis- 
tribution in terms of the shape parameters 

_ _ _ v „oo 2(l+2j) ^ nOC 3(2j + l) 

+YC . (3+AQ-D) ( E1 °) 

Finally, the normal stress balance (|E9[) yields c jm o. 

Expanding around sphere, for the linear flow Q 
Eq. Q yields 



dt 



- c JOT 2 + /jm at r = 1 



(Ell) 



The first term represents the motion of the interface due 
to the normal component of the velocity and the second 
term describes the rotation of the deformed shape. 

Substituting Cj m 2 in Eq. (|E11[) yields the evolution 
equation ([32]) with the following coefficients 



C(X,j,m) = d(XJ)- 1 [cf m V7(7 + T) (2j + 1) 
+cf m2 (4j 3 + 6j 2 -4 ? -3)] 



and 



where 



and 



r(A,cr ,i) = -E(j,a ) 



d(A,j) 



d(A,j)-9A- 1 + (-5 + 3j 2 + 2j 3 ) 



E(j» = (j + 2) (j - 1) (i(j + 1) + a) 



(E12a) 
(E12b) 

(E12c) 
(E12d) 



For a simple shear flow (J33J), solving Eq. (j3"2")) for the sta- 
tionary amplitudes in the tank-treading regime gives 



/2±2 = T" 



i X 10, 



E(2,a )Tix 
and for the inclination angle 



(9+23A) 
6 



1 



arctan 



X(9 + 23A) 



(E13) 



(E14) 



These expressions agree with Seifert's results for a vesicle 
with no viscosity contrast (A = 2) 14 . At long times, the 
excess area is stored in the /2±2 modes only. Substituting 
the tension jF4| in (|E12dp leads to 



(23A + 9) 



<-l 



1920tt 



(23A + 9) 2 A' 



(E15) 



Thus the dependence of the mode amplitudes on the cap- 
illary number in Eq. (|E13|) cancels. Expanding Eq. (|E14|) 
for small values of the excess area A we obtain (|4U|) . 

The surface solenoidal velocity field u^ ml satisfies iden- 
tically the incompressibility condition (fT6)) . The stress 
boundary condition 



r° m \- (A- l)r™ 1= 



(E16) 



corresponds to a spherical viscous drop with no surface 
tension, i.e. the u^ nl flow field is unaffected by the in- 
terfacial stresses. Consequently, for the j = 1 component 
of the shear flow, which corresponds to rigid body rota- 
tion, we have that C\ m \ = c^ nl , i.e. at leading order the 
particle rotates with the flow. 

For the sake of completeness we mention the relation 
between our notation and the one used by Seifert 1 ^. 



Cjm2 A-j'm 7 

CjmO = {Yjm + 2Xj m ) + 1)]~ 



(E17) 
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APPENDIX F: THE ISOTROPIC TENSION 



and 



The area constraint (|28|) serves to determine the 
isotropic part of the tension ct - We can split Eq. (|E12b|) 
into 



/3(A,i) = -2a(j) 



j(j + 1) 



(F3) 



where 



r(A,cr ,i) = +(T /3(\,j) . 



ot{\j) = -2a(j)- 



(Fl) 



d(X,j) ' 



The modes have to satisfy the area constraint A = 0. 
Multiplying the evolution equations (j32jl by a(j)f* m , 
(F2) summing up and solving for the tension leads to 



c 



Ejro a 0') [C( A ' J> m )fjm + X 1 a{X, j)f jm f. 



X 1 Ejm a (j)/ 3 ( A I j)/jm/ 



(F4) 



The complicated dependence of the tension on the shape 
modes makes the evolution equations Eq. (|32p nonlinear. 
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